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Abstract. Atomic and molecular data for the transitions of a number of astrophysically interesting species 
are summarized, including energy levels, statistical weights, Einstein ^-coefficients and collisional rate coef- 
ficients. Available collisional data from quantum chemical calculations and experiments are extrapolated to 
higher energies (up to E/k ~ 1000 K). These data, which are made publically available through the WWW 
at http : //www . strw . leidenuniv . nl/ ^moldata are essential input for non-LTE line radiative transfer programs. 
An online version of a computer program for performing statistical equilibrium calculations is also made available 
as part of the database. Comparisons of calculated emission lines using different sets of collisional rate coefficients 
are presented. This database should form an important tool in analyzing observations from current and future 
(sub) millimetre and infrared telescopes. 
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A wide variety of molecules has been detected in space 
to date ranging from simple molecules like CO to 
more complex organic molecules like ethers and alcohols. 
Observations of molecular lines at millimetre and infrared 
wavelengths, supplemented by careful and detailed mod- 
elling, are a powerful tool to investigate the physi cal and 
chem i cal condition s of astrophysical objects (e.g.. iGenzell 
Il99ll iBlackl lioOO). To constrain these conditions, lines 
with a large range of critical densities and excitation tem- 
peratures are needed, since densities typically range from 
~ 10 2 — 10 9 cm~ 3 and temperatures from ~ 10 — 1000 K 
in the interstellar and circumstellar environments probed 
by current and future instrumentation. 

In recent years, different molecules have been devel- 
oped as tracers for different physical and chemical con- 
ditions (see van Dishoeck & Hogerheijde 1999 for a re- 
view). For example, CO is used as a tracer of the total 
gas mass whereas readily observed molecules with large 
dipole moments, such as CS, HCO + and HCN constrain 
the density structure. The wide variety of H2CO and 
CH3OH lines accessible at millimetre and submillimetre 
wavelengths trace both the temperature and density struc- 



ture (e.g.. lMangum fe Wootten1ll993(l . Organic molecules 
like CH3OCH3 and CH3CN probe the chemical complex- 
ity. Deuterated molecules contain a record of the condi- 
tions and duration of the cold pre-stellar phase. Si- and 
S-bearing molecules, in particular SiO and SO2, probe 
shocks. Lines of the main species as well as the (gener- 
ally) optically thin isotopes are needed to determine ac- 
curate abundances and line profiles. High frequency lines 
and vibrationally excited lines are particularly valuable for 
probing the warm and dense inner parts of the circumstel- 
lar en velopes (e.g.. lZiurvs fe Turnerll98rilBoonman et alJ 
l200lh . 

To extract astrophysical parameters, the excitation 
and radiative transfer of the lines need to be calcu- 
lated. Indeed, it is becoming increasingly clear that more 
information - including chemical gradients through- 
out the source — can be inferred from the data if 
a good molecular excitation mode l is available (e.g., 
ISchoier et all2002l iMaret et all2004j) . The simplest mod- 
els adopt the 'local' approximation, for example in the 
widely used large velocity gradient (LVG) method. A 
number of more sophisticated, non-local radiative trans- 
fer codes have been develope d for the int e rpretation of 
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a review). The application of these codes ranges from pro- 
tostellar environments to the circumstellar envelopes of 
late-type stars. The radiative transfer analysis requires ac- 
curate molecular data in the form of energy levels, statisti- 
cal weights and transition frequencies as well as the spon- 
taneous emission probabi lities and collisiona l rate coeffi- 
cients. Th e JPL 1 catalog jPickett et alJll99sh . HITRAN 2 
datab ase faothman et al.ll2003h . and the CDMS 3 cata- 
logue l|Muller et al.ll200ljr contain energy levels and tran- 
sition strengths for a large number of molecular species. 
Detailed summaries of the theoretical methods and the 
uncertainties involved in determining collisio n al rat e co- 
efficients are given by iGreenl l)l975a|) . iRouefll (^990) and 
iFlowerl l|l99fj|) . In this paper, these and other literature 
data on the rotational transitions of 23 different molecules 
are summarized and extrapolations of collisional rate coef- 
ficients to higher energy levels and temperatures are made. 
The molecular data files can be found at the webpage 
http://www.strw.leidenuniv.nl/~moldata and is the 
first effort to systematically collect and present the data 
in a form easily used in radiative transfer modelling of 
interstellar regions. The focus is on rotational transitions 
within the ground vibrational state, but the lowest vibra- 
tional levels are included for a few common species where 
such data are availabl e. Many of the data files presented 
here were adopted by ISchoier et alJ l)2002f) to model the 
circumstellar environment of the protostar IRAS 16293- 
2422. In addition, data files for three atomic species are 
presented. The excitation of atomic fine structure levels 
plays an important role in cooling of a wide variety of 
astrophysical objetcs. 

An online version of RADEX 4 , a statistical equilib- 
rium radiative transfer code using an escape probabil- 
ity formalism, is made available for public use as part of 
the database. RADEX is comparable to the LVG method 
and provides a useful tool for rapidly analysing a large 
set of observational data providing constraints on phys- 
ical conditions, such as density a nd kinetic temperature 
l|jansen et alJ 11994 Ijansenlll995|) . RADEX provides an 
alternative to the w idely used rotatio n temperature di- 
agram method (e.g., iBlake et al.lll987h which relies upon 
the availability of many optically thin emission lines and is 
useful only in roughly constraining the excitation temper- 
ature in addition to the column density. A guide for using 
the code in practice is provided at the RADEX homepage. 
RADEX will be presented in more detail in a forthcoming 
paper (van der Tak et al., in prep.) at which point the 
source code will be made publically available. 

2. Energy levels 

In this section the molecular structure is briefly reviewed. 
This serves merely to provide some basic information 

1 http://spec.jpl.nasa.gov 

2 http : //cf a-www . harvard . edu/HITRAN/ 

3 http://www.cdms.de 

4 http : //www . strw . leidemmiv . nl/ ~moldata/radex . html 



needed to properly use the data files. Detailed discussions 
on molecular (and atomic) structure can be found in, e.g., 
iTownes fc Schawlowl lll975l) . 

2.1. General considerations 

The energy levels are obtained from the JPL, HITRAN, 
and CDMS catalogues. The energy levels and the corre- 
sponding line frequencies are thus of spectroscopic quality 
and may be used for the purpose of line identification, 
unless stated otherwise. 

Generally, we retain only the ground vibrational 
state and include energy levels up to E/k ~ 1000 K. 
Vibrationally excited levels are usually not well popu- 
lated in the regions probed by current (sub)millimetre 
telescopes. Moreover, little is known about colli- 
sional rate coefficients fo r vibrational transitions (e.g., 
IChandra fc Sharmal l200ll) . However, for some specific 
molecules, e.g. CO and CS, vibrationally excited levels 
are also included. In the prolongation of this project, data 
files including vibrational levels will be added for more 
molecular species. 

Molecules with ortho and para versions (or A- and E- 
type as in the case of e.g. CH 3 OH) are treated as separate 
species. 

2.2. Linear molecules 

The energy levels for diatomic and linear polyatomic 
molecules in the 1 E electronic state are quantified, to first 
order, according to 

E = BJ(J+l), (1) 

where B is the rotational constant and related to the mo- 
ment of inertia /, around axes perpendicular to the inter- 
nuclear axis, through B = (2/) _1 . Heavy linear molecules, 
like HC3N, have more densely spaced energy levels than 
diatomic molecules like, e.g., CO. These pure rotational 
energy levels are classified according to the rotational 
quantum number J and their statistical weights are 

S=(2J + 1). (2) 

Note that to obtain state energies of spectroscopic accu- 
racy, Eq. Q must be augmented with centrifugal distor- 
tion (tx J 2 (J + l) 2 ) and higher-order terms. The majority 
of molecular species presented here have a S electronic 
ground state, i.e., the sum of the orbital angular momenta 
of their electrons and the sum of the electron spins are 
both zero. However, there are some exceptions where ei- 
ther can be non-zero. 

For a molecule in a 2 £ electronic ground state, e.g., SO 
and CN, the sum of the electron spins is 1/2. The non-zero 
spin creates a splitting of the levels due to coupling be- 
tween the electron spin and the total angular momentum 
of the molecule. The total angular momentum is quantified 
according to N and includes the rotation of the molecule. 
Molecules like, e.g., O2 have a total electron spin of 1 in 
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a 3 S electronic ground state. Some important molecules 
such as NO, NS, and OH have a 2 H ground state with a 
total electronic orbital momentum of 1 and total spin of 
1/2. Spectroscopically, such molecules show 'A— doubling', 
with 2 Ui/ 2 and 2 n 3 / 2 ladders. 

The various molecular angular momenta may couple 
together in many different ways, such as spin-orbit and 
spin-spin coupling. Ideally, these fall in one of five differ- 
ent classes, known as Hund's coupli ng cases. In practice, 
interm ediate cases often occur; see iTownes fc Schawlowl 
GjffiJ for details. 



2.3. Non-linear molecules 

The structure of non-linear molecules, such as e.g. H 2 CO 
and CH3OH, is more complex. Rotation can take place 
around different axes of inertia, characterized by the ro- 
tational constants A, B and C which, in absence of any 
symmetry, involve different amounts of energy, A>B>C . 
The degree of asymmetry is measured by Ray's parameter 

2B-A-C 
K= A-C ' (3) 

and is —1 for a prolate symmetric top (B=C, e.g. CH3CN) 
and +1 for an oblate symmetric top (B=A, e.g. NH3). 
Asymmetric rotors such as H 2 have \k\ <C 1. 

The energy levels of symmetric top molecules, such as 
NH3 and CH3CN, are described by the quantum numbers 
J and K, where K is the projection of the total angular 
momentum J on the symmetry axis. For a prolate sym- 
metric top molecule, the energy of a rotational level is 
given (to first order) by 



E = BJ(J + 1) + (A - B)K 2 



(4) 



The energy levels for a slightly asymmetric prolate top 
such as H 2 CO can be calculated from 



B + C 



J(J + 1) + [A 



B + C 



(5) 



where w~K 2 with correction s due to the slight asymmetry 
Townes fc Schawlowl l|l975h . Appendix III]. 



2.4. Hyperfine splitting 

A further complication arises when the nuclear spin cou- 
ples to the rotation producing what is known as hyper- 
fine splitting. The astrophysically most relevant cases are 
when the molecule contains a 14 N or D nucleus. When 
the lines are spectroscopically resolved, hyperfine struc- 
ture provides information on the optical depths, which i s 
otherwise hard to obtain fe.g.. lSchmid-Burgk et alJfefJO^ . 

Hyperfine splitting can be important in line transfer 
and introduce non-local effects for lines overlapping 
in frequency, e.g., the J = 1 — » lines of N2H 4 " and 
HCN. The common assumption is that the hyperfine 
components each have the same excitation temper- 
ature. However, exceptions to this rule ('hyperfine 



anomalies') have been observed and more d e tailed - 
treatments developed llStutzki fc Winnewisserl Il985t 
| Truong-Bach fc Nguven-Q-Rieul Il989t ILindavist et alJ 
2000). Often the splitting between individual hyperfine 
components is small, producing lines which are separated 
in frequency by a small amount compared to the line- 
broadening, so that this splitting can be safely neglected 
and treated as a single level for the purpose of excitation 
analysis. 

The first release of the database includes hyperfine 
splitting for some of the most relevant molecules, such 
as HCN and OH. Future releases will present data files 
with hyperfine splitting included for additional species. 

3. Radiative rates 

3.1. General formulae 

The radiative rates for dipole transitions from an upper 
state u to a lower state I can be calculated from 



An, 



64tt 4 z/V S ul 



3c 3 h 



fju 



(6) 



where /x is the electric dipole moment and S is the tran- 
sition strength. The transition strength depends on the 
complexity of the molecule and is explained below in some 
detail. Strictly speaking, the dipole moment should be av- 
eraged over the vibrational wavefunction(s) of the tran- 
sition involved (fJ- v ), but in practice often the dipole mo- 
ment appropriate for the equilibrium geometry is taken 
(He)- The electric dipole moments are assumed to be the 
same for all isotopes of a particular molecule, even though 
small differences exist for fj, v . Because of the v 3 factor, the 
resulting Einstein A-values can still differ considerably for 
isotopes, especially for deuterated species. 

For transitions with AJ±1 in linear molecules, the 
transition strength is 



(7) 



whereas for symmetric top molecules, the transition 
strength from level J, K to J — 1, K is given by 



S u i = 



J 2 -K 2 
J 



(8) 



In the general case of asymmetric tops, simple expressions 
for S do not exist. 



3.2. Dipole moments 

Transition strengths are available from the spectroscopic 
databases mentioned above (JPL, CDMS, HITRAN). 
There is some inconsistency in the astrophysical litera- 
ture regarding the choice of values of electric dipole mo- 
ments, however. This often manifests itself as an appar- 
ent bias against results of ab initio theoretical calcula- 
tions, even when experimental results for transient species 
are merely estimated or wholly absent. A case in point 
concerns the pair of ions HCO + and HOC + : the widely 
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cited JPL catalogue offers [i = 3.30 D for HCO + and 
fi = 4.0 D for HOC" 1 " based o n low -level theoretical es- 
timate s of IWoods et al ] dl975l) and iGudeman fc Woodsl 
ill 9821). respectively, wherea s accurate ab initio values 
from lBotschwina et alJ l|l993|) gi ve fx n (HCO+) = 3 9 3±.01 
D and /x e (HOC+) = 2.74 D. IZiurvs fc Annonil ijl 995ft 
adopted a similar val ue, Ato(HOC + ) = 2.8 D, from an ab 
initio computation of iDefrees et alJ (|l982l) . Because the 
inferred column densities scale as oc [i~ 2 , these discrepan- 
cies in dipole moments can result in errors of factors of 
two in derived abundances. 

Table n collects values of dipole moments for a (non- 
exhaustive) sample of molecules of astrophysical interest. 
Users are encouraged to remain aware of the original lit- 
erature. Unless otherwise indicated, all entries refer to the 
electronic and vibrational ground states. 

For small dipoles, centrifugal corrections to the dipole 
moment are appreciable. In the case of CO, rotational 
effects reduce the A-value by 1% for J=7 and by 10% for 
J=22. The JPL and CDMS catalogues consider this effect 
and so do our datafiles. 

4. Collisional rate coefficients 

4.1. General considerations 

The rate of collision is equal to 

Cui = n co i7 ui , (9) 

where n co \ is the number density of the collision part- 
ner and "f u i is the downward collisional rate coefficient (in 
cm 3 s _1 ). The rate coefficient is the Maxwellian average 
of the collision cross section, er, 

{8kT^ ~ 1/2 



Table 1. Summary of adopted dipole moments"'' 



lul 



1 

k~T 



^ aEe- E / kT dE, 



(10) 



where k is the Boltzmann constant, \x is the reduced mass 
of the system, and E is the center-of-mass collision energy. 
The upward rates are obtained through detailed balance 



Jlu = lul — e ' 
9i 



(11) 



where g is a statistical weight. 

The collisional rate coefficients -f u i usually pose the 
largest source of uncertainty of the molecular data input 
to the radiative transfer analysis (however, see discussion 
on dipole moments in Sect. 13.2]) . The dominant collision 
partner is often H2 except in photon dominated regions 
(PDRs) where collisions with electrons and H can become 
important. The collisional rate coefficients presented here 
are mainly with H2 and only in a few cases (in partic- 
ular the atoms) are collisions with H and electrons also 
treated. Where available, the data files include collisions 
with ortho- and para-H2, e.g., in the case of CO. 

If only data for collisions with He are available, a first 
order correction can be made by assuming H2 to have 
the same cross sections. This approximation is strictly 



Molecule 3 
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Ue 


Method 


Reference 




[D] 


[D] 






CO 


0.110 




— expi 


I 


SO 


1 52 + 02 




GXpt . 


2 




1.633 




expt. 


3 




1 Q5R 4- 005 




exp 


4 


SiO 


3.098 
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SiS 


1.73 ± .06 




expt' 


g 


HCO + 


— 1 cm + 01 


— S 90 + 01 


UiU hi hhhbYJ 


7 


HOC + 


2.74 


2.8 


nn iTiifid 

Lb w hi LhiihYJ 


8 9 


OPS 


n 71 59 




expt. 


i n 


HC3N 


3.732 






xx 


HCN 


2.985 




expt' 


12 


HNC 


3.05 ±.1 




expt. 


13 


c-C 3 H 2 


3.27 ± .01 




expt. 


1 


CH 3 CN 


3.922 




expt. 


14 


H 2 CO 


2.332 




expt. 


15 


N 2 H+ 




3.4 


ab initio 


16 


HCS+ 


1.958 




ab initio 


17 


CH3OH 


0.896 




expt. 


18 


NH 3 


1.476 ± .002 




expt. 


19 


H 2 


1.847 




expt. 


20 


SiC 2 


2.393 ± .006 




expt. 


21 


HC1 


1.109 




expt. 


22 


OH 


1.655 




expt. 


23 


H 3 0+ 


1.44 




ab initio 


24 



a Same data are also adopted for isotopes and deuterated 
species, unless stated in the datafile. 

b All values are in units of debye (D), where 1 D = 10 -18 esu 
cm. When the original source has presented /j, in atomic units, 
a conversion fa ctor of 1 au = eap = 2.5417 5 D has been 
Refs - m iGoorvitchl Jl994h. (2) lLovas et al. 
(3) iPatel et all Jl979l) . (4) IWinnewisser fc Cook 
(5) [Ravmonda et al] \l97(t) '. (6) [Hoeft et alf 
(7) [Botschwina et alJ (ll993^ . ( 8) iBotschwinal 
(91 Defrees et alJ Jl98aT (10) iMuenteil il96gT). 



1992) 
1968) 
196c) 
1989) 



iDeLeon fc Muenterl 



,.(12) 





Ebcnstcin & Mucntcr 



dl984). (131 iBlackman et all Jl976t) ; (14) iGadhi et alJ jll995h 
(15) iFabricant et all lll97^ 



iBotschwina &i Sebald ) 
(19) ICohen fc Povnteil 

1 J1985I 



C amv- Pevre^^^J 




(16) iGreen et all il974l). (17 ) 
(18) ISastrvetalJ dl98lD . 

(20) From JPL based on 



De Leeuw fc Dvmanua , 
(24) IBotschwina et al.l(fl98 



(21) ISuenram et~ai] (I198 1 
lll97lTl . (23) iPeterson et alJ 



only valid for very cold sources, where most H2 is in the 
ground J — state without angular momentum. Then 
from Eq. (|10|> the rate coefficient for collisions between a 
molecular species X and H 2 

/ \ 1/2 

7x-h, = 7x-Hc — — - I ■ (12) 

\MX-H 2 / 

If the mass of the molecule is much larger than that of He 
and H 2 , the scaling factor is 1.4. 

Some molecules of significant interest lack calculated 
collisional rate coefficients. In these cases the rates for a 
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similar molecule have been adopted and only scaled for 
the difference in reduced mass following Eq. (|12|l . This 
procedure works best for O— >S substitutions (for exam- 
ple, scaling HCO + rates for the case of HCS + ) since such 
molecules have a similar molecular structure. 

For most species, only rate coefficients with He or H2 
J—0 are available. Values with H2 J=l can be larger by 
factors of 2-5 due to supplementary terms in the interac- 
tion potential (e.g. Green 1977, H 2 example). This addi- 
tional uncertainty is often not considered in astrophysical 
analyses. In the case of CO and H 2 0, separate rate coeffi- 
cients are available for collisions with ortho- and para-H2. 
The online version of RADEX weighs these coefficients 
by the thermal value of the H2 o/p-ratio at the kinetic 
temperature. The o/p-ratio is approximated as the J=l 
to J=0 population ratio with a maximum of 3.0, which 
is an overestimate by at most 20% (at T—155 K). In the 
datahles available for download the collisional rate coef- 
ficients for collisions with ortho-H2 and para-H2 are kept 
separate. 

To obtain the collision rate, RADEX simply multiplies 
the collisional rate coefficients with the H2 density. To in- 
clude the effect of collisions with He, the user must multi- 
ply the density by 1.14 (to first order) for a He abundance 
with respect to H2 of 20%. 

The adopted collisional rate coefficients are presented 
in Tables [21 and [3] for atomic and molecular species, re- 
spectively. For isotopomers the same set of collisional rate 
coefficients as for the main isotope was adopted, unless 
otherwise stated. Tables and |3 show the temperature 
range and maximum energy (-E max ) for which calculations 
are available. Also, the collision partner is indicated. Only 
the downward values arc given in the data files; the upward 
rate coefficients are obtained through detailed balance us- 
ing Eq. dHJ). 



4.2. Accuracy of adopted rate coefficients 

Most of the collisional data summarized in Table 2 have 
been obtained from theoretical calculations, with exper- 
imental cross-checks possible for only a few cases. Most 
experiments reflect the average of many collisional events, 
with comparisons typically done for relaxation rates and 
collision-induced pressure line broadening. State-to-state 
measurements have been possible for only a few systems 
and they often report relative rather than absolute cross 
sections. Also, experiments with H2 are usually done for n- 
H 2 (i.e., 75% o-H 2 and 25% p-H 2 ), rather than for H 2 J=0 
or 1. Nevertheless, such comparisons between theory and 
experiment, as well as those between different theoretical 
methods, have given some indication of the uncertainties 
in the collisional rate coefficients. Excellent accounts of the 
methods involved and details on ind ividu al systems ar e 
given hv lCreenl ill 975ah . lFklwel lIlOTflh andlRouefll |l9P,r| : 
recent developments are reviewed bv lRoueff et all l|2004j) . 
Here only a brief summary is given. 



Table 2. Summary of atomic collisional data from the 
literature. 



Atom 




T 
[K] 


-E'max 

[cm" 1 ] 


Collision 
partner 


Ref. 


C 


10 


- 200 


43 


H 


1 




10 


- 20000 


43 


e _ 


2 




100 


- 2000 


43 


H+ 


3 




10 


- 150 


43 


He 


4 




10 


- 1200 


43 


H 2 


5 


C+ 


5 


- 3162 


63 


H 


1 




10 


- 20000 


43054 


c _ 


6 




10 


- 250 


63 


H 2 


7 





50 


- 1000 


227 


H 


1 




50 


- 3000 


227 


c _ 


8 




100 


- 100 


227 


H+ 


9 




20 


- 1500 


227 


H 2 


10 



Refs. - m lLaunav fc Rouefj il977f> . f2l Ijohnson et alJ jl987Tl . 

(31 iRoueff fc Le Bourlol Jl99Ch . (41 IStaemmler fc Flowed 
il99lTl. (5) ISchroder et all il991f> . (6 1 I Wilson fc Belli faOoK 
(71 iFlower fc Launavl Jl977|). (81 iBell et all Il998l> . (91 



IChambaud et alJ Jl98oh . (lOl ljaauet et all (119921 



The theoretical determination of collisional rate coeffi- 
cients consists of two steps: (i) determination of the inter- 
action potential V between the colliding systems; and (ii) 
calculation of the collision dynamics. Significant progress 
in the second part has been made in the last decades, 
aided by the increased computer speed. The most accu- 
rate method is the Close-Coupling (CC) method, in which 
the scattering wave function is expanded into a set of ba- 
sis functions. This method is exact if an infinite number 
of basis functions or 'channels' is taken into account. In 
practice a finite number of channels is used, resulting in 
a set of coupled second-order differential equations. The 
absolute accuracy of the results can easily be checked by 
increasing the basis set, and is of order a few % for a given 
interaction potential. This method works very well for low 
collision energies and relatively light species, although care 
should be taken at the lowes t energies whether resonance s 
are properly sampled (e.g.. iDubernet k, Grosieanl 12002^1 . 
However, the method becomes increasingly computation- 
ally demanding at high energies and for heavy polyatomic 
molecules with small splittings between the rotational en- 
ergy levels resulting in many channels. 

The most popular approximate dynamical methods 
are the Coupled States (CSt) or 'centrifugal decoupling' 
method and the Infinite Order Sudden (IOS) approxima- 
tion. In the CSt method, the centrifugal potential is as- 
sumed to conserve the projection of the angular momen- 
tum on the axis perpendicular to the plane of the collision 
partners. This approximation is often valid at higher en- 
ergies if the collision is dominated by the repulsive part 
of the potential. In the IOS approximation, an additional 
assumption is that the molecule does not rotate during 
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Table 3. Summary of molecular collisional data from the 
literature and new extrapolated rate coefficients. 
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collisions. This may be appropriate for heavy rotors at 
energies much larger than the rotational energies. From 
comparisons with the more exact CC results, it is found 
that absolute uncertainties for the CSt method range from 
^10% to a factor of 2, with lesser uncertainties in the rela- 
tive values. The propensities in the collisions are recovered 
correctly. In contrast, the IOS method can have uncer- 
tainties up to an order of magnitude. Computer programs 
which include the CC, CSt and IOS options are publically 
available (see Hutson & Green 1994 5 , Flower et al. 2000 6 , 
Manolopoulos 1986 and Alexander & Manolopoulos 1987 7 

The above quoted ranges of uncertainties assume that 
the interaction potential is perfectly known. Often, this is 
not the case and the potential surfaces form the largest 
source of error in the collisional rates with uncertain- 
ties that are difficult to assess. The interaction potential 
consists of a short-range repulsive part, an intermediate- 
range interaction part where a weak molecular bond is 
formed, and a long-range part dominated by electrostatic 
interaction. The intermediate part is most difficult to de- 
termine and requires high-level quantum chemical mod- 
els. The most accurate method is that of Configuration 
Interaction (CI), but it can become very costly in com- 
puter time. Other methods include Hartree-Fock Self- 
Consistent-Field (SCF) and perturbation methods, and 
more recently Density Functional Theory (DFT), but each 
of these methods has its drawbacks. An old approximate 
method, the Electron Gas model, is now obsolete, but 
some dynamics calculations for astrophysical systems still 
use these potentials (e.g., CS-H 2 , Turner et al. 1992). 

The following selected examples serve to illustrate the 
range of absolute errors in the adopted collisional rate 
coefficients. It should be noted that relative values often 
have less uncertainty and that these are most relevant for 
astrophysical applications: small absolute errors can often 
be compensated by small adjustments in the abundance 
of the species. 

4.2.1. CO-H 2 




Early calcul at ions by IGreen Thadd^ul (1 
ISchinke et all l|l985iT and iFlower fc Launavl |l 

illustrate the sensitivity of the results to different poten- 
tial energy surfaces. Absolute differences in individual 
collisional rate coefficients range from a few % up to 
40%, with the relative values usually having less scatter. 
Comparison of com puted cross sections u s ing a new 
CO-H 2 potential bv Ijankowski fc Szalewic^ l)l998(l with 
pressure broaden ing and scattering experiments by 
iMengel et alJ l)2000|) suggests an overall average absolute 
accuracy of better than 10% at T > 30K, but somewhat 
less good at the lowest temperatures where the deviations 
can increase to 30-50%. No information is available on 
the accuracy of the larger A J transitions (e.g., A J > 10), 



http : / / www . giss . nasa . gov/molscat 

:tti1 



http : / / ccp7 . dur . ac .uk/molcol .html 
7 http : //www. chem.umd. edu/physical/ alexander/hibridon 
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Fig. 1. Predicted CO line intensities, using different sets 
of calculated collisional rate coefficients, for an isothermal 
homogeneous sphere with a kinetic temperature 10 K, a 
H2 density of 10 3 cm -3 and a CO column density of 3 x 
10 16 cm~ 2 .The line intensities are shown in relation to the 
values o btained using the CO-pH 2 rate coefficients from 
iFlowerl (|200l4 - The upper rotational quantum number 
J u is indicated on the s-axis. The rotational transitions 
are out of thermal equilibrium and, for transitions below 
J = 4 — > 3, optically thick. 



which become important at high temperatures such as 
found in dense shocks. The same potential surface has 
been u s ed in t he latest set of rate coefficients given by 
IFlowerl lj2001al) which are adopted here. 

The following simple test problem illustrates the con- 
sequences of using different sets of collisional rate coef- 
ficients. Line intensities were calculated for the lowest 5 
rotational transitions of CO for a molecular cloud of con- 
stant temperature and density using RADEX. The model 
has a temperature of 10 K, H2 density of 1 x 10 3 cm -3 
and a total CO column density of 3 x 10 16 cm -2 over a 
line width (full- width at half-maximum) of 1 km s _1 . All 
lines are out of thermal equilibrium and the three lowest 
rotational transitions are optically thick. As is shown in 
Fig. n differences of up to ±150% are found, especially 
for collisions with para-H2 compared with ortho-H 2 . 



4.2.2. H2CO-H2 

The H2CO-H 9 rate coeffici ents given in our database are 
obtained from Grcc n1 lll991^ who calculated values for the 
Hs CO-He system using a very old potential energy surface 
bv lCarrison et alJ ( 1975^ based on SCF and limited CI cal- 
culations. These rate coefficients and the adopted surface 
have recently been tested against pressure broadening and 
time-resolv ed double-resonance studie s for three low-lying 
transitions ijMengel fe De Lucial2000|) . Satisfactory agree- 
ment is found for the H 2 CO-He system, with differences in 
cross sections ranging from a few % up to 20%. The devia- 
tions are largest at the lowest temperatures, <10 K, as was 
also found for CO-H2. For H2 as the collision partner, the 
cross sections are found to be up to a factor of two higher, 



significantly more than the value of 1.4 expected from the 
difference in masses, illustrating that simple scaling from 
He collisions may introduce errors up to 50%. 

4.2.3. OH-H2 

One of the computationally most challenging systems is 
OH-H2, since OH is an open shell molecule with a 2 n 
ground state so that two potential surfaces and hyperfinc 
splitting are involved. Resul ts for collisions with both o-H? 
and P-H2 are presented by Offer fe van DishoeclJ l|l992f) 
using an old potential surf ace base d on SCF calculat ions 
l|Kochanski fc Flowerll98lL and bv lOffer et all lll994h us- 
ing a new surface computed using CT (|Offer fc van Hemerti 
Il993h . The differences due to the potential energy sur- 
face range from 10 % to more than an order of mag- 
nitude for individual rate coefficients. Comparison with 
state-to-state experimental cross sections with both n-H2 
and P-H2 at one specific energy gives surprisingly good 
agreement, usually within 50% but with occasional excur- 
sions up to an order of magnitude 1 Schreel fc ter Meulenl 
1996(). Moreover, all the propensities for individual hyper- 
fine transitions are well reproduced. 



4.3. Adopted collisional rate coefficients 

Below follows a summary of the collisional rate coefficients 
adopted in the first release of the database. Molecules for 
which only one set of calculated collisional rate coefficients 
is available and where no extrapolation was performed 
are not described further here. The principle method for 
extrapolating the downward collisional rate coefficients 
(AJ=J u ^Ji, J u >Ji) in temperature in the case of a lin- 
ear molecule is l)de Jong et al.lll975l iBieging et al.l ll998) 



7uJ =a(AJ)yexp[-b(AJ)y L ^] x exp[-c(AJ)y L ^], (13) 

where y = AE u i/kT and the three parameters a, b, and 
c are determined by least-squares fits to the initial set 
of rate coefficients for each AJ . This reproduces most of 
the calculated rate coefficients to within 50% and typically 
within 20%. For more details on the extrapolation scheme, 
including extrapolation in energy levels, see Sect. [5] (only 
available in the online version of this journal). 

When the specified kinetic temperature falls outside 
the region where collisional rate coefficients are available, 
i.e. from T\ ow to Thigh, RADEX makes no further extrapo- 
lation and assumes the downward rate coefficients at Ti ow 
and Thigh, respectively. 

4.3.1. CO 

For CO the c ollisional rate coefficients calculated by 
IFlowerl l 2001ah have been adopted as a starting point. 
These computations cover temperatures in the range from 
5 K up to 400 K and include rotational levels up to J = 29 
and J = 20 for collisions with para-H 2 and ortho-H 2 , re- 



8 



F. L. Schoier et al.: An atomic and molecular database for analysis of submillimetre line observations 



■ Flower (2001) ^=1^0 
-a Flower & Launay (1985) 

• Schinke et al- (1985) 

- + Larsson et al. (2002) 
& This work 

- □ This work 




100 1000 





10 100 
T [K] 





100 1000 



i<r 12 










itr 13 




/=40a0n 


icr 14 






i<r 16 






lO" 16 






lO" 17 






io- la 






io- 19 






io- 20 






IO" 21 






10 22 






io- 23 






io- 24 






IO" 25 






, -26 







100 1000 



Fig. 2. Calculated and extrapolated collisional de-excitation rate coefficie nts for CO in collis i ons w ith para-tb. Open 
triangles indicate extrapolation in temperature to the rate coefficients of iFlower fc Launavl l)l985j) (filled triangles). 
O pen squares sho w the extrapolation to higher temperatures and energy levels of t he recent rate coefficie nts calculated 
bv lFlowerl l)2001a|) (filled squares). F or comparison the rat e coefficients presented bv lSchinke et~al] lll985h (filled circles) 
and the extrapolation performed bv lLarsson et alJ l)2002|) (filled stars) are shown. 



spectively. Both sets of rate coefficients were then extrapo- 
lated to include energy levels up to J = 40 [using Eq. I|18|) ] 
and temperatures up to 2000 K [using Eq. (|13[) ] , as de- 
scribed in Sect. 16. ll fn the datafile available for download, 
the collisional rate coefficients for collisions with ortho-H.2 
and para-H2 are kept separate. However, in RADEX they 
are weighted together as described in Sect. 14. ll 

Fig. El shows the extrapolation of CO collisional de- 
excitation rate coefficients for collisions with para-El^. It 
is clear that extrapolated rate coefficients are uncertain 
and depend on both the original data set from which the 
extrapolation is made and the method adopted. However, 
the extrapolated values typically agree within 50% in the 
case of CO. The largest discrepancies, up to an order of 
magnitude, naturally arise in the region where extrapola- 
tion in both temperature and energy levels are performed. 
Thus, in the parts of parameter space where extrapo- 
lated rates are being used to infer physical conditions, care 
should be taken as to any astrophysical conclusions drawn 
from the modeling. 



in Sect. 16. ll No extrapolation to temperatures lower than 
20 K was attempted. 

4.3.3. SiO 

For S iO the rate coefficients calculated by iTurner et alJ 
l|l992|) have been adopted, computed for temperatures in 
the range 20 — 300 K and including rotational levels up 
to J = 20 for collisions with H2. This set was then ex- 
trapolated to include energy levels up to J = 50 [using 
Eq. (HHJ)] and temperatures up to 2000 K [using Eq. (fT3|) ]. 
as described in Sect. 16. ll No extrapolation to temperatures 
lower than 20 K was attempted. 

4.3.4. SiS 

No calculated rate coefficients are available for SiS. 
Instead, the same set of collisional rate coefficients as for 
SiO has been adopted. 



4.3.2. CS 



For C S the rate coefficients calculated by ITurner et alJ 
(1992) have been adopted as a starting point. These val- 
ues have been computed for temperatures in the range 
20 — 300 K and include rotational levels up to J = 20 
for collisions with H2. This set was then extrapolated to 
include energy levels up to J = 40 [using Eq. l(T%|) ] and 
temperatures up to 2000 K [using Eq. (|T3"|) ]. as described 



4.3.5. HCO+ 

The rate coefficien ts for HCO + in collisions with H2 have 
been calculated bv iFlowerl 1 19991) for temperatures in the 
range 10 — 400 K and rotational levels up to J = 20. This 
set of rate coefficients was then extrapolated to include en- 
ergy levels up to J = 30 [using Eq. (|T%|) ] and temperatures 
up to 2000 K [using Eq. Q], as described in Sect.lO 
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4.3.6. HC 3 N 

The rate coefficien ts for HC.^N in colli s ions w ith He have 
been calculated bv lGreen fc Chapman! ljl97S|) for temper- 
atures in the range 10 — 80 K and rotational levels up to 
J = 20. This set of rate coefficients was then extrapolated 
to include energy levels up to J = 50 [using Eq. lfT%|) ] and 
temperatures up to 2000 K [using Eq. 1|13|)]. as described 
in Scct. lbMl The rate coefficients were then scaled by 1.39 
to represent collisions with H 2 instead of He. 

4.3.7. HCN 

The rate coefficien ts for HCN in collisions w ith He have 
been calculated bv lGreen fc Thaddeusl l)l974j) for temper- 
atures in the range 5 — 100 K and rotational levels up 
to J = 7. This work has subsequently been extended 
by S. Green (unpublished data) to include rotational lev- 
els up to J = 29 and temperatures from 100 — 1200 K. 
Extrapolation of the rate coefficients to include energy 
levels up to J = 29 for temperatures below 100 K [us- 
ing Eq. Ijl8(l ]. as described in Sect. 16.11 has been made. 
The rate coefficients were subsequently scaled by 1.37 to 
represent collisions with H 2 instead of He. The collisional 
rate coefficient s between various hy p erfine levels have been 
calculated by iMonteiro fc Stutzkl ljl986[) for the lowest 
(J < 4) rotational levels and temperatures from 10 — 30 K 
in collisions with He. A datafile based on these collisional 
rate coefficients is also made available separately. 

4.3.8. HNC 

No calculated rate coefficients are available for HNC. 
Instead, the same set of collisional rate coefficients as for 
HCN has been adopted. 

4.3.9. N 2 H+ 

The rate coefficients for N2H" 1 " in collisio ns with He atoms 
have been calculated bv iGreenl l|l975bj) for temperatures 
in the range 5 — 40 K and rotational levels up to J = 6. 
Given the limited range in temperature and energy levels, 
we have instead adopted the same rate coefficie nts as for 
HCO + . This is motivated by the discussion in IMonteirol 
l)l984|) where the rate coefficients for these two species in 
collisions with He are found to be very similar, typically 
within 10%. 

4.3.10. HCS+ 

The rate coefficients for HCS + in c ollisio ns with He atoms 
have been calculated bv IMonteirol l|l984|) for temperatures 
in the range 10 — 60 K and rotational levels up to J = 
10. This set of rate coefficients was then extrapolated to 
include energy levels up to J = 23 [using Eq. l(T%|l ] and 
temperatures up to 1000 K [using Eq. (|13|) ]. as described 
in Sect. 16. ll The rate coefficients were subsequently scaled 
by 1.38 to represent collisions with H 2 instead of He. 



4.3.11. H 2 

In RADEX the r ate coefficients for H2O in collisions with 
He calculated bv lGreen et alJ l|l993|) are used as default. 
The rates were computed for temperatures in the range 
from 20 to 2000 K including energy levels up to about 
1400 cm -1 . These rate coefficients were subsequently 
scaled by 1.35 to represent collisions with H 2 instead of 
He. In addition, a datafile containing the recent rate co- 
efficie nts for H?Q in collisions with P-H? llGrosiean et all 
120031) and o-H 2 (|Dubernet fc GrosiearJ l2002h calculated 
for low temperatures (5 — 20 K) has been constructed. In 
the datafiles available for download, the rate coefficients 
for collisions with ortho-H 2 and para-H 2 are kept sepa- 
rate. However, in RADEX they are weighted together as 
described in Sect. 5.1. 

4.3.12. S0 2 

For non-linear molecules there are no simple scaling re- 
lations such as Eq. I|13(l . In Sect. 6.2 (only available in 
the online version of this paper) the procedure adopted 
to extrapolate rate coefficients for S0 2 is presented. As 
starting point the calculated ra te coefficients for S0 2 in 
collisions with He calculated bv iGreerJ l|l995|) were used. 
These rates were computed for temperatures in the range 
from 25 to 125 K including energy levels up to about 
62 cm" 1 . Extrapolation was made to include energy levels 
up to 250 cm -1 and temperatures in the range from 10 to 
375 K. The rate coefficients were subsequently scaled by 
1.4 to represent collisions with H 2 instead of He. 

5. Summary 

A compilation of atomic and molecular data in a homo- 
geneous format relevant for radiative transfer modelling 
is presented. The data files are made available through 
the WWW and include energy levels, statistical weights, 
Einstein A-coefficients and collisional rate coefficients. 
Extrapolation of collisional rate coefficients are generally 
needed and different schemes for this are reviewed. 

In addition to the atomic and molecular database, an 
online version of a computer code for performing statis- 
tical equilibrium calculations is made available for use 
through the WWW. The program, named RADEX, is an 
alternative to the widely used rotation diagram method 
and has the advantage of supplying the user with physical 
parameters such as density and temperature. 

Databases such as these depend heavily on the efforts 
by the chemical physics community to provide the relevant 
atomic and molecular data. We strongly encourage further 
efforts in this direction, so that the current extrapolations 
of collisional rate coefficients can be replaced by actual 
calculations in future releases. 
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6. Extrapolation of collisional rate coefficients 

6.1. Linear molecules 

An often adopted starting point when fitting and extrap- 
olating collisional rate coefficients is to take advantage of 
the IOS approximation in which the entire matrix of state- 
to-state rate coefficients can be calculated from the basic 
7lo rate coefficients (e.g. iGoldflam et al"lll977j) 



J+J' 

L=\J-J'\ 



1) 



J J' L 




where 



J J 1 L 




7lo, (14) 



(15) 



is the Wigner 3-j symbol. This expression is valid only 
in the limit where the kinetic energy of the colliding 
molecules is large compared to the energy splitting of the 
rotational levels. Since the energy splitting increases with 
J this expre ssion becomes less acc urate for higher rota- 
tional levels. IPePristo et alJ l|l979f) show that by multi- 
plying Eq. 11411 (within the summation) with 

A(X,J)- 6 + ^ 2 



where 



O(J') = 0.13 J'B I 



1/2 



(16) 



(17) 



one can approximately correct for this deficiency. Here 
Bq is the rotational constant in cm" 1 , I is the scattering 
length in A (typically I w 3 A), [i is the reduced mass of 
the system in amu and T is the kinetic temperature in K. 
Extrapolation of the rate coefficients down to the lowest 
J=0 level can be made both in temperature as well as 
in J allowing t he general state-to-state coefficient s to be 
extended (e.g.. lAlbrecblll983l:rLa"rsson et al.ll2002^ . 

Alternatively, and in line with the IOS approximation, 
the downward collisional rate coefficients (AJ=J u ^Ji, 
Ju>Ji) can be extrapolated in temperature using Eq. I|13l) . 
Given its simplicity we have adopted this procedure for 
extrapolation of the rate coefficients in temperature. 

Extrapolation to include higher rotational levels was 
carried out by fitting the collisional rate coefficients con- 
necting to the ground rotational state, at a particular tem- 
perature, to a second order polynomial 

7j = exp(a + bJ + c,J 2 ), (18) 

where a, b and c are parameters determined from the fit. 
Figure |3] illustrates the fit to collisional rate coefficients 
down to the ground rotational state for CO-H2 using 
Eq. Ijl8|) . Similar extrapolations can be made in temper- 
ature. However, here we have adopted t he ap proach by 
Ide Jong et alJ l)l975|) and iBieging et all ()l998|) and used 
Eq. (|13f) for the extrapolation in temperature. This ex- 
tends the fit over a larger range of energies. The IOS 
approximation [Eq. (|14fl ] was then used to calculate the 
entire matrix of state-to-state rate coefficients. The CO 
molecule is used in Sect. l4.3T1 to illustrate the above men- 
tioned schemes. 
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Fig. 3. The solid line s are fit s to the CO-P-H2 collisional 
rate coefficients from iFlowerl l|2001a|) (squares ) for tran- 
sitions down to the ground state from upper energy levels 
J u using a second order polynomial. 



6.2. Non-linear molecules 

For non-linear species there are no simple scaling relations 
and one has to resort to custom-made fitting formulae for 
each ca se. The only case conside red here is that of SO2 
used bv Ivan der Tak et al J l)2003h . In the prolongation of 
this project extrapolated collisional rate coefficients will 
be presented for additional non-linear molecules. As the 
starting p oint fo r inela stic collisional data for SO2, the 
results of iGreenl (|l995() were used. However, those data 
only cover the lowest 50 states, up to 62 cm" 1 (J w 12), 
while states up to J = 25 are commonly observed. 

Figure 0] plots Green's downward rate coefficients, 
summed over all final states, as functions of initial state. 
These sums approach asymptotic values for E u > 40 cm" 1 . 
Deviations from this behaviour due to detailed quantum 
mechanical selection rules are seen not to exceed 20%. The 
figure also shows that the rate coefficients increase approx- 
imately as T 1 / 2 , again to «20% accuracy. This behaviour 
indicates that the total rate coefficients only depend on 
temperature through the collision velocity, while the de- 
excitation cross sections are constant. 

Fig. 03 shows that most collisions lead to de-excitation 
into states that are not far down in energy. The thick black 
curve is our fit to this behaviour: it is the normalized mean 
of the various thin light curves which represent Green's 
data. Transitions by more than 15 states are considered 
negligible. 

Based on these trends, the rate coefficients for de- 
excitation of SO2 in inelastic collisions with He are ex- 
trapolated as follows. For the 50 lowest states, Green's 
values at 25 < T < 125 K are used and multiplied 
by (T/125K) 1 / 2 at temperatures up to 375K and down 
to 10 K. For states between 62 and 250 cm -1 above 
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Fig. 4. Collisional de-excitation rate coefficients for the 
lowest 50 states of SC^summed over all lower levels, cal- 
culated from the data by Green (1995) for various tem- 
peratures. 



ground, a total de-excitation rate coefficient of 1.0 x 
cm -3 s -i j g assumec i ; shown by Fig. 0]to be a 
good zeroth-order description for other levels. The state- 
to-state rate coefficients are derived by multiplyin g these 
totals by the mean propensities from lGreenl 1 1995|) . given 
by the black curve in Fig. [SJ All results are multiplied 
by 1.4 to account for the mass difference between H2 and 
He. While this procedure is admittedly crude and does 
not take the detailed quantum mechanics of the interac- 
tion into account, it catches the spirit of more detailed 
calculations. 



Fig. 5. State-to-state de-excitation rate coefficients for 
SO2 as fractions of the total downward rate coefficient 
(Fig. , as a function of the number of levels by which the 
transition is changed. The light (coloured) curves are val- 
ues from Green (1995) at T=25 K for the 10 th , 20 th , 30 th , 
40 th and 50 th state above ground. The thick black curve 
is the normalized mean of the light (coloured) curves, 
adopted here to extrapolate Green's rate coefficients to 
higher- lying levels. The states are labelled in order of in- 
creasing energy. 



